QCD Sum Rules for the X(3872) as a mixed molecule-charmoniun state 
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We use QCD sum rules to test the nature of the meson A (3872), assumed to be a mixture between 
charmonium and exotic molecular [cg][gc] states with J PC — 1 ++ . We find that there is only a small 
range for the values of the mixing angle, 0, that can provide simultaneously good agreement with 
the experimental value of the mass and the decay width, and this range is 5° < 6 < 13°. In 
this range we get m x = (3.77 ± 0.18) GeV and T(X -> J/im+Tv') = (9.3 ± 6.9) MeV, which 
are compatible, within the errors, with the experimental values. We, therefore, conclude that the 
X(3872) is approximately 97% a charmonium state with 3% admixture of ~88% D°D*° molecule 
and ~12% D + D*~ molecule. 

PACS numbers: 11.55.Hx, 12.38.Lg , 12.39.-x 



I. INTRODUCTION 

Among the new hadronic states discovered in the last 
few years, the AT (3872) is one of the most interest- 
ing. It has been first observed by the Belle collabora- 
tion in the decay B+ -> X(3872)K+ -> J/^tt+tt" K+ 
This observation was later confirmed by CDF, DO 
and BaBar 2]. The current world average mass is 
m x = (3871.4 ± 0.6) MeV which is at the threshold 
for the production of the charmed meson pair D°D°*. 
This state is extremely narrow, with a width smaller 
than 2.3 MeV at 90% confiedence level. Both Belle and 
Babar collaborations reported the radiative decay mode 
AT (3872) -> 7 J/> @,[3, which determines C = +. Fur- 
ther studies from Belle and CDF that combine angular 
information and kinematic properties of the tt + tt~ pair, 
strongly favor the quantum numbers J PC = 1 ++ or 2 

I SI- 

In constituent quark models [7| the masses of the pos- 
sible charmonium states with J PC — 1 ++ quantum num- 
bers are: 2 3 Pi (3990) and 3 3 Pi (4290), which are much 
bigger than the observed mass. In view of this large 
mass discrepancy the attempts to understand the X me- 
son as a conventional quark-antiquark states were aban- 
doned. The next possibility explored was to treat this 
state as a multiquark state, composed by c, c and a 
light quark antiquark pair. Another experimental find- 
ing in favor of this conjecture is the the fact that the de- 
cay rates of the processes X(3872) — > J/ip 7r + 7r _ 7r° and 
X(3872) — > J/-07r + 7r _ are comparable [3j: 



X -> J/V'TT+TT-TT 

X -> J/ifjn+n- 



= 1.0 ±0.4 ±0.3. 



(1) 



This ratio indicates a strong isospin and G parity vi- 
olation, which is incompatible with a cc structure for 
X(3872). 
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In a multiquark approach we can avoid the isospin vi- 
olation problem. The next natural question is: is the 
X made by four quarks in a bag or by a meson-meson 
molecule? 

The observation of the above mentioned decays, plus 
the coincidence between the X mass and the D*°D° 
threshold: M(D*°D°) = (3871.81 ± 0.36) MeV 0, in- 
spired the proposal that the A (3872) could be a molec- 
ular (D*°D° - D*°D°) bound state with small binding 
energy @, [3. The D*°D° molecule is not an isospin 
eigenstate and the rate in Eq. ^ could be explained in a 
very natural way in this model. 

Maiani and collaborators [llj suggested that A (3872) 
is a tetraquark. They have considered diquark- 
antidiquark states with J PC = 1 ++ and symmetric spin 
distribution: 

x q = [cq]s=i[cq]s=o + [cq]s=o[cq)s=i- (2) 
The isospin states with 1 = 0, 1 are given by: 

X u + X,i Wt ... A„ — X d 



X(I = 0) 



V2 



X(I=1) 



V2 



(3) 
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In [ll| the authors argue that the physical states are 
closer to mass eigenstates and are no longer isospin eigen- 
states. The most general states are then: 

Xi = cos6X u +s\n6X d , X h — cos9X d -sm8X u , (4) 

and both can decay into 2-7T and 3n. Imposing the rate in 
Eq.Q, they obtain 6 ~ 20°. They also argue that if A/ 
dominates B + decays, then Xh dominates the B° decays 
and vice-versa. Therefore, the X particle in B + and B° 
decays would be different with [HEI M(X h )-M(Xi) = 
(8±3) MeV. There are indeed reports from Belle [13| and 
Babar *14j Collaborations on the observation of the B° — > 
K° X decay. However, these reports (not completely 
consistent with each other) point to a mass difference 
much smaller than the predicited ~ 8 MeV. 

All the conclusions in ref. were obtained in the 
context of a quark model. Given the uncertainties inher- 
ent to hadron spectroscopy, it is interesting to confront 
these theoretical results with QCD sum rules (QCDSR) 
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calculations. This was partly done in [15] where, us- 
ing the same tetraquark structure proposed in ref. [ill ], 
the mass difference M{Xh) — M(Xi) was computed and 
found to be in agreement with the BaBar measurement 
(M(X h ) - M(Xi) = (3.3 ± 0.7) MeV). The same cal- 
culation [lH has obtained mx — (3.92 ± 0.13) GcV. In 
QCDSR we can also use a current with the features of the 
mesonic molecule of the type (D*°D° - D*°D°). With 
such a current the calculation reported in [l6| obtained 
the mass mx = (3.87 ± 0.07) GeV in a better agree- 
ment with the experimental mass. Therefore, from a 
QCDSR point of view, the A (3872) seems to be bet- 
ter described with a D*D molecular current than with 
a diquark-antidiquark current. We feel though that the 
subject deserves further investigation. 

In this work we use again the QCDSR approach to 
the X structure including a new possibility: the mixing 
between two and four-quark states. This will be imple- 
mented folowing the prescription suggested in [171 ] for the 
light sector. The mixing is done at the level of the cur- 
rents and will be extended to the charm sector. In a 
different context (not in QCDSR), a similar mixing was 
suggested already some time ago by Suzuki [1^|. Phys- 
ically, this corresponds to a fluctuation of the cc state 
where a gluon is emitted and subsequently splits into a 
light quark- antiquark pair, which lives for some time and 
behaves like a molecule-like state. As it will be seen, in 
order to be consistent with X decay data, we must con- 
sider a second mixing between: (cc) + (D*° D° — D*° D°) 
and (cc) + {D*+D- - D*~D + ). 

With all these ingredients we perform a calculation of 
the mass of the A (3872) and its decay width into 27r and 
3tt. 



II. THE MIXED TWO-QUARK / FOUR QUARK 
OPERATOR 

There are some experimental data on the A (3872) me- 
son that seem to indicate the existence of a cc component 
in its structure. In ref. [HI it was shown that, because 
of the very loose binding of the molecule, the production 
rates of a pure A (3872) molecule should be at least one 
order of magnitude smaller than what is seen experimen- 
tally. Also, the recent observation, reported by BaBar 
[H, of the decay A(3872) -> ip(2S)-y at a rate: 



ff(A^(2ff) 7 ) 
B(X -> iry) 



3.4 ± 1.4, 



is much bigger than the molecular prediction [191 ]: 

Mil) ~ 4 x to-. 
T(X - V>7) 



(5) 



(6) 



While this difference could be interpreted as a strong 
point against the molecular model and as a point in favor 
of a conventional charmonium interpretation, it can also 
be interpreted as an indication that there is a significant 



mixing of the cc component with the D°D*° molecule. 
Similar conclusion was also reached in refs. [2(| [H|. 
Therefore, we will follow ref. [l7| and consider a mixed 
charmonium- molecular current to study the A (3872) in 
the QCD Sum Rule framework. 

For the charmonium part we use the conventional axial 
current: 



(7) 



The D° D*° molecule is interpolated by 22, 23, 24]: 



(x) = 4= 
V2 



(u a (x)^ 5 C a {x)cb(x)j fl U b (x)) 



(u a (x)j fJ ,C a (x)c b (x)j 5 U b (x)) 



(8) 



As in ref. [17| we define the normalized two-quark current 
as 



A2u) = 



1 



6V2 



(9) 



and from these two currents we build the following mixed 
charmonium-molecular current for the A(3872): 

j;(x) = Bin(9)j^(x) + cos(0)jf<>(*). (10) 



III. THE TWO POINT CORRELATOR 



The QCD sum rules [2J, |27J, |28| are constructed from 
the two-point correlation function 



U MV {q) = 1 / d A x e*- x (0|T[J?( a: )J?t(o)]|o) 



Q 2 



n (<z 2 ) 



q 2 ' 



(11) 



As the axial vector current is not conserved, the two func- 
tions, IT! and n , appearing in Eq. (fTTj) are independent 
and have respectively the quantum numbers of the spin 
1 and mesons. 

The sum rules approach is based on the principle of 
duality. It consists in the assumption that the correla- 
tion function may be described at both quark and hadron 
levels. At the hadronic level (the phenomenological side) 
the correlation function is calculated introducing hadron 
characteristics such as masses and coupling constants. At 
the quark level, the correlation function is written in in 
terms of quark and gluon fields and a Wilson's opera- 
tor product expansion (OPE) is used to deal with the 
complex structure of the QCD vacuum. 

The phenomenological side is treated by first 
parametrizing the coupling of the axial vector meson 
1 ++ , X, to the current, J", in Eq. (fT0|) in terms of the 
meson-current coupling parameter X u : 



\j;\x) 



X u e L 



(12) 
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Then, by inserting intermediate states for the meson X, 
we can write the phenomenological side of Eq. (TTTj) as 



(A 



u\2 



•>X 



-5V 



9m 1v 



x 



(13) 



where the Lorentz structure projects out the 1 ++ state. 
The dots denote higher mass axial-vector resonances. 
This ressonances will be dealt with through the intro- 
duction of a continuum threshold parameter sq- 

In ref. [2!| it was argued that a single pole ansatz can 
be problematic in the case of a multiquark state, and 
that the two-hadron reducible (2HR) contribution (or S- 
wave DD* contribution, in the present case) should also 
be considered in the phenomenological side. However, in 
ref. [3(| it was shown that the 2HR contribution is very 
small. The reason for this is the following. The 2HR 
contribution, in our case, can be written as 13011: 



d 4 p ( -gjq> +p tl y v /m 2 1 



(2tt)4 



(p - i) 2 - m 2 D 



where 



(0\j^\DD*(p)) = \ DD »s fl (p). 



(14) 



(15) 



Following ref. 30] the current two-meson coupling: 
Add*, can be written in terms of the D meson decay 
constant, /d, and the coupling of the D* meson with a 
4-quark current. This last quantity should be very small, 
because the properties of the D* meson, both in spec- 
troscopy and in scattering, are very well understood if 
it is an ordinary quark-antiquark state. Therefore, the 
parameter Add* , should be very small, as in the case of 
the pentaquark [30(, and the 2HR contribution can be 
safely neglected. 

In the OPE side we work up to dimension 8 at the 
leading order in a s . The light quark propagators are cal- 
culated in coordinate-space and then Fourier transformed 
to the momentum space. The charm quark part is cal- 
culated directly into the momentum space, with finite 
m c , and combined with the light part. The correlator in 
Eq. (fTTj) can be written as: 



iW?) 



\6V2 



cos 



in(26»)) n$*\q) 



6V2 



(16) 



with: 



4V } (?) 



(fx e* 



\T\j®(x)jM (Q)]\Q). (17) 



After making a Borel transform of both sides, and 
transferring the continuum contribution to the OPE side, 



the sum rule for the axial vector meson up to dimension- 
eight condensates can be written as: 

(A U )V 

+ M(sin ( 2^))nr)(M 2 ) + 



+ sm 2 (0)n[ 4A \M 2 ) 7 



(18) 



where: 



n^(M 2 )=/ ds e-^%% 2 J t (s) +U%%(An (19) 



llll 2 



n^\M 2 ) = [° ds e-^ P ^( s ) +ngl >( M*), 

J4m? 



(20) 



U[ 4 - 4 \M 2 ) = / dse- s ' M2 

J 4m 2 



Ppert\ S ) + P(uu)\- 



[(44) 

(uu) (uGu) 



.( 44 ) fo\ _L ^(44) 



+ 



+ n£*> r am 2 ), 

(uu) (uGu) V / ' 



and 



(22) , s s ( 4m 2 ^ 2 



(21) 



(22) 



nf G 2 >(M 2 ) 



__(1&X rl d 
3-2 5 tt 2 JO u " 



2a(l-q)M 2 +TO 2 
a(l~a)M 2 



+2m 



2 (2a-l)m 2 c +a(a 2 -1)M 2 
c M 4 a 3 (a-1) 



e °<i-°)m 2 , 



o (24) (s) - 



^o (22) ( S ) 

Hpert\°/> 



6V2 



(23) 



(24) 



rf 24 ) __Mtt( 22 ) 
6V2 



n^>(^) = -fin^ ) (M^) (25) 



U (2i l AM 2 ) = 5{q9 ;% G f fidaj£-e-«l=5z*, (26) 



l-(q+/3) 2 y 4 , 



(27) 
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P(«4 (s) = J amin daj 0min d(3 y -p >-K 2 (a,(3), 14 

(28) 



P(G 2 ) \ I -TTT^fB-J " a J/3 mi „ 



4m 2 



(29) 



9 (l-(a+/3) ) 



l-2a-2/3 



(30) 




o (44) fsl 



3m e 

" 2 S 7T 4 



(it<7<7 • Gm) f" ma: * da 



2(m^-a(l-a)s) 
1-a 



(31) 



n| 44) w ^ ,(m 2 ) 



m^(un) (ug s <7-Gu) rl 



2 5 7T 2 



Jo da 



a(l-a)M 2 +ro 2 
a(l-a)M 2 
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FIG. 1: Relative contributions of the terms in eqs. (1221) to 
(|32[) grouped by condensate dimensions. We start with the 
perturbative contribution and each subsequent line represents 
the addition of one extra condensate dimension in the expan- 
sion. 
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and we define K(a, f3) = (a + j3)m 2 — a(3q 2 . 

By taking the derivative of Eq. (fT8"|) with respect to 
1/M 2 and dividing the result by Eq. (fT8|) we can obtain 
the mass of mx without worrying about the value of 
the meson-current coupling X u . The expression thus ob- 
tained is analised numerically using the following values 
for quark masses and QCD condensates [HI, HH : 

m c (m c ) = (1.23 ±0.05) GeV, 
(uu) - (0.23±0.03) 3 GeV 3 , 

4( 



(uga.Gu) = m^uu), 

ml = 0.8 GeV 2 , 
(g 2 G 2 ) = 0.88 GeV 4 . 



(33) 



In Fig. [T] we show the contributions of the terms in 
Eqs. (j2"2")l to (|32p grouped by condensate dimensions di- 



vided by the RHS of Eq. JT8|). We have used s 



1/2 



4.4 



GeV and = 9°, but the situation does not change much 



Ivr(GeV') 

FIG. 2: The dashed line shows the relative pole contribution 
(the pole contribution divided by the total, pole plus con- 
tinuum, contribution) and the solid line shows the relative 
continuum contribution. 



for other choices of these parameters. It is clear that the 
OPE is converging for values of M 2 > 2.6 GeV 2 and we 
will limit our analysis to that region. 

The upper limit to the value of M 2 comes by imposing 
that the QCD pole contribution should be bigger than the 
continuun contribution. The maximum value of M 2 that 
satisfies this condition depends on the value of so, being 
more restrictive for smaller sq. In Fig. [5] we show a com- 
parison between the pole and continuun contributions 

1/2 

for the smaller s$ we will be considering (s = 4.4) and 
9 = 9°. The condition obtained from Figgis M 2 < 3.2 
GeV 2 , but in this case, the dependence on the choice of 
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9 is very strong. Taking into account the variation of 9 
we have determined that, for 5° < 9 < 13°, the QCDSR 
are valid in the following region: 



2.6 GeV 2 < M < 3.0 GeV^ 



(34) 



In Fig. [3l we show the X meson mass in this region. 
We see that the results are reasonably stable as a function 
of M 2 . From Fig. [3] we obtain m x = (3.80 ± 0.08) GeV 




Nf(GeV') 

FIG. 3: The X meson mass as a function of the sum rule 
parameter (M 2 ) in the region of eq. (|34[) for different values of 



the continuum threshold: sj^ 2 = 4.4 GeV (solid line), Sq /z = 
4.5 GeV (dashed line) and s\ /2 = 4.6 GeV (dotted line). 

where the error includes the variation of both s$ and M 2 . 
If we also take into account the variation of 9 in the region 
5° < 9 < 13° we get: 



1/2 



m x = (3.77 ±0.18) GeV, 



(35) 



which is in a good agreement with the experimental 
value. The value obtained for the mass grows with the 
value of the mixing angle 9, but for 9 > 30° it reaches a 
stable value being completely determined by the molec- 
ular part of the current. 

From Eq. (fT5)) we can also obtain A" by fixing mx 
equal to the experimental value (mx = 3.87 GeV). Using 
the same region in 9, sq and M 2 that we have used in 
the mass analysis we obtain: 



A" = (3.6±0.9).10~ d GeV . 



(36) 



IV. 



DECAY OF THE X(3872) AND THE THREE 
POINT CORRELATOR 



As discussed in Sec. I, one of the most intriguing facts 
about the meson A" (3872) is the observation, reported 
by the BELLE collaboration [3], that the X decays into 
J/ip 7r + 7T~7r°, with a strength that is compatible to that 



of the J/-07T+7T - mode, as given by Eq .|T]). This de- 
cay suggests an appreciable transition rate to J/ip to and 
establishes strong isospin violating effects. It still does 
not completely exclude a cc interpretation for X since 
the origin of the isospin and G parity non-conservation 
in Eq. ([I]) could be of dynamical origin due to p° — u> mix- 
ing [32]]. However, the observation of the ratio in Eq. ([I]) 
is an important point in favor of the molecular picture 
proposed by Swanson [l9[ . In this molecular picture the 
A(3872) is mainly a D°D* a molecule with a small but 
important admixture of pJ/i/j and wJ ' components. 

It is important to notice that, although a D°D*° 
molecule is not an isospin eingenstate, the ratio in Eq. ([1]) 
can not be reproduced by a pure D°D*° molecule. This 
can be seen through the observation that the decay 
width for the decay X — > J/tpV — > J/ipF where F = 
7r + 7r~ (7r + 7T~7r°) for V = p(u>) is given by [Til. 1331] 



as 



■W) 



1 



8nm x 



\M\ 2 B V ^ F 



p(s) 



(s-m 2 ,) 2 + {m v T v f 



where 



p(s) = 



yj\{m x ,ml,s) 



2m x 



(37) 



(38) 



with A(a, 6, c) = a 2 + b 2 + c 2 - lab - 2ac - 26c. The 
invariant amplitude squared is given by: 



l-^l 2 = 9x^vf{ m x,m^, 



(39) 



where gx^v is the coupling constant in the vertex 
XJ/ipV and 



f{m x ,m^,s) = - Am x h - 



2s 



(m 2 x -s) 2 \ m 2 x -m 2 , + s 



2m% 



2m\ 



Therefore, the ratio in Eq. (JTJ) is given by: 

T(X -> J/ipir+Tr~ir°) _ gx^m^Y ^B^^l^ 



T(X -> Jf-ipTT+n ) g 2 Xipp m p T p B p ^^I p 
where 



(40) 



(41) 



Iv = I dslf(m x ,m,p,s) 

x P(s) 

(s - m 2 v ) 2 + (m v TvY 



(42) 



Using B^,„ = 0.89, T w = 8.49 GeV, m w 
782.6 MeV, B p ^„ = 1, T p = 149.4 GeV and m p 
775.5 MeV we get 



T(X -> J/ijj 



Y{X J/iP 



£pQ = 0.118 (^V. (43) 

b 7T*7T ) \9XtPpJ 



6 



The couplings, gx^v-, can be evaluated through a 
QCDSR calculation for the vertex, X(3872) J/ipV, that 
centers in the three-point function given by 

IWp,p',<z) = J <?xd*y e ip '- x e^m^fay), (44) 
with 

lWz, y ) = {0\T\j+(x)tf{y)g\0)]\0), (45) 

where p = p' + q and the interpolating fields are given 
by: 



■ ip _ - 



(46) 



Ju = ^-{UaluUa + {-l) Iv d alv d a ), (47) 



with N p = l,_Ip = 1, N„ = 1/3 and I„ = 0. If X(3872) 
is a pure D°D*° molecule, is given by Eq. (jSJ). In this 
case the only difference in the OPE side of the sum rule 
is the factor Ny and, therefore, regardless the approxi- 
mations made in the OPE side and the number of terms 
considered in the sum rule one has 

n£„ a (p, j/, q) = N v Il°„ p a E (p,p', q). (48) 

To evaluate the phenomenological side of the sum rule 
we insert, in Eq. (|45|) . intermediate states for X, J/ip and 
V. We get [13: 



n(fjrW><?) = 



iXxm^f^myfy gx^v 
(p 2 - m x )(p' 2 - m 2 , )(q 2 - m\) 



( - e a " va {p' a + q„) - ^nts^L 
V m v 



m l )' 



(49) 



Therefore, for a given structure the sum rule is given 



by: 

i\xm^,f^m v f v gx^v 
(p 2 - m 2 x )(p' 2 - m^iq 2 - m 2 v ) 



N v Il OPE (p,p',qK>0) 



from where, considering m p ~ one gets: 

f u _ N u _ 1 
gxippfp N p 3 

Using f p = 157 MeV and f u = 46 MeV we obtain 



gx^ 



1.14, 



and using this result in Eq. (|4"3"| we finally get 

T(X -► J/^TT+TT-TT ) 



T(X -> J/^TT+TT-) 



0.15. 



(51) 



(52) 



(53) 



It is very important to notice that this is a very gen- 
eral result that does not depend on any approximation 
in the QCDSR. This result shows that the admixture of 
pj /ij) and ujJ/ip components in the molecular model of 
ref. [1 91 ] is indeed very important to reproduce the data in 
Eq. |T]). It is also important to notice that, in a QCDSR 
calculation of the decay rate X — > J/ipV, the cc admix- 
ture in the D°D*° molecule, as given by Eq. (|10[) . does 
not solve the problem of geting the ratio in Eq. |T]) . This 
can be seen by using, in Eq. (|45j) . = J™, with J™ given 
by Eq. flTUJ). One gets: 



sin(0)n™°<(x,y), 



(54) 



where 



K^y) = {oiTVtWfitof^ioW), (55) 



and 



n^(x,y) = (OlT^^^jl^^WllO), (56) 

with j„ 2 ^ and ji 4 "' given by Eqs. J7]) and ©. Using the 
currents in Eqs. (|4"Tj) and (|4"6"| for the mesons U and J/tp, 
it is easy to see that 



Kl a (x,y) 



N v 



Tr hp,S c ac (x) lal5 S c ca (-x)} x 



x Tr [ 7 „S&(0) + (-1)^7,^(0)] .(57) 

For V = p with J p = 1 the result in Eq. (|57|) is obvi- 
ously zero due to isospin conservation, in the case that 
the quark u and d are degenerate. However, even for 
V = u) (I u = 0), the result in Eq. ([57]) is zero be- 
cause Tr [7 M ,Sf 6 (0)] = 0. Therefore, in the OPE side, 
the three-point function is given only by the molecular 
part of the current in Eq (fT0|) : 



n fiVa (x,y) = S ia(9)Rl n ^(x,y) 



(58) 



that can not reproduce the experimental observation in 
Eq. HJ, as demonstrated above. 

In the following, to be able to reproduce the data in 
Eq.([T]), instead of the admixture of pJ/i/> and cuJ/ip com- 
ponents to the D°D*° molecule, as done by Swanson 
[ijj, we will consider a small admixture of D + D*~ and 
D~D* + components. In this case, instead of Eq. fTU)) we 
have 



= cos a J" (x) + sinaJ^(x), 



(59) 



with J£(x) and J°(x) given by Ea.([l0")). 

If we consider the quarks u and d to be degenerate, 
i.e., m u — rrid and (uu) = (dd), the change in Eq. (|10p 
to Eq. ([59|) does not make any difference in the results in 
Sec. III. 
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By inserting , given by Eq. (|59|) . in Eq. (|45f and 
considering the quarks u and d to be degenerate, one has 

^-ij, V a{p,p',q) = sin((9)^=(cosa 

+ (-l)^mna)ll^(p,p', q ), (60) 

with 

n°„ p a E (p,p', q ) = J d+uj ( Tr lVMOtsx 

x Sl bl {-y) lv Sl lb {y) la S c bal {k~p')\ + 
- Tr [l,S c a , c (k) la Sl b , {-y)l„St, b {y) lb S c ba , (k - 



In the phenomenological side, considering the defini- 
tion of A" in Eq. (fT2"|) and the definition of the current in 
(f59|) . we can define 

Xx = cos a\ u + sin a\ d — (cos a + sin a)A q , (62) 

where X q was evaluated in Sec. Ill, and is given in 
Eq. (j3"6"|) . Using Eq. (j6"2"| in Eq. P5|) , the phenomenological 
side of the sum rule is now given by: 



m l )' 



(63) 



From Eqs. (|60|) and (|63|) we get the following relation 
between the coupling constants: 



gxjnju _ ^(cosa + sin a) 
gxippfp N p ( cos a — sin a) 



(64) 



Using the previous result in Eq. (|4"3"]) and the numerical 
values for f w and / p we have 



T(X -> J/^TT+TT-TT ) 
r(X J/lpTT+TT-) 



0.15 



cos a + sm a 



cos a — sm a 



(65) 



i(cosa + sin a)\ q m^f^mvfv gx-ipv 
(p 2 - m 2 x )(p /2 - m^jiq 2 - m 2 v ) 



x -e 



This is exactly the same relation obtained in refs. [Ill.l33l]. 
that determines a ~ 20° for reproducing the experimen- 
tal result in Eq.(|T}. A similar relation was obtained in 
ref. [HI where the decay of the X into two and three 
pions goes through a D D* loop. 

With this mixing angle a defined, we can now eval- 
uate the decay rate itself, for any one of the decays: 
X — » J/ipp or X — > J/tpu), since they will be the same. 
Therefore, we choose to work with X — > J/tpoj since the 
combination cos a + sin a appears in both sides of the 
sum rule and the result for gx^u is independent of a. 




FIG. 4: Diagrams which contribute to the OPE side of the sum rule. 



In the OPE side we consider condensates up to di- 
mension five , as shown in Fig. [¥] Taking the limit 
p2 _ p/ _ pi anc j doing a single Borel transform to 
P 2 — > M 2 , we get in the structure e^^p'^q^p'^ (the same 



considered in ref.fH) (Q 2 = -q 2 )- 

C{Q 2 ) (e- m V M2 - e" m i/ M2 ) + B e~ a »l M * = 

(Q 2 + m 2 )n(°^)(M 2 ,Q 2 ), (66) 
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where 



shown bellow: 



U (OPE \M 2 ,Q 2 ) = 



6V2tt 2 Q 2 
du e~ u/M2 ^l-Am 2 / 



3Q 2 " 
1 m 2 



Ami 



ci = 2.5 x 10~ 2 GeV 7 , 

c 2 = 38 GeV 2 , 
B = 2.9 x 1(T 4 GeV 5 . 



(70) 



TO n f 1 , 1 + 3a 2lL_ 

— / da e«( 1 -°) J,! 

16 J a 



In Eq. (J66 



C(Q 2 ) 



sin(0) 



rriufu 



(67) 



:9x^{Q 2 ), (68) 



In Fig.[5]we can see that the Q 2 dependence of C(Q 2 ) 
is well reproduced by the chosen parametrization in the 
interval 2.5 < Q 2 < 4.5 GeV 2 , where the QCDSR is 
valid. 



and B gives the contribution of the pole-continuum tran- 
sitions |33l . l35l . [36j . so ancl u o are the continuum thresh- 
olds for X and J/-0 respectively. Notice that in Eq. ([67|) 
we have introduced the form factor gx^{Q 2 )- This is 
because the meson uj is off-shell in the vertex XJ/tpco. 





FIG. 



Q 2 (GeV 2 ) 

Momentum dependence of C(Q 2 ) 



for s 



1/2 



: 4.4 GeV. 



The solid line gives the parametrization of the QCDSR results 
(dots) through Eq. JBS} and (JTOJ). 

The form factor gxipu{Q 2 ) can then be easily obtained 
by using Eqs. (f6"5| and (f6"9"|) . Since the coupling constant 
is defined as the value of the form factor at the meson 



pole: Q 2 = 



to determine the coupling constant we 



have to extrapolate gx^(Q 2 ) to a Q 2 region where the 
sum rules are no longer valid (since the QCDSR results 
are valid in the deep Euclidian region). Using = 
3.1 GeV, ra x = 3.87 GeV, fa = 0.405 GeV, A" = 3.6 x 
FIG. 5: Values of C(Q 2 ) obtained by varying both Q 2 and 10 3 GeV 5 from Eq. J3H]) and varying in the range 5° < 



3.5 3.0 
M 2 (GeV 2 ) 



M in Eq. ((66. 



If we parametrize C(Q 2 ) as a monopole: 



9 < 13°, we get: 



gxi>u = gxipuji-ml) = 5.4 ± 2.4 



(71) 



C{Q 2 ) = 



ci 



Q 2 + c 2 ' 



(69) 



The decay width is given by: 

r {X -» J/^+TT-TT ) = 



we can fit the left hand side of Eq. (|6"o]) as a function 
of Q 2 and M 2 to the QCDSR results in the right hand 
side, obtaining c\, ci and B. In Fig. [5] we show the 
points obtained if we isolate C(Q 2 ) in Eq. (|66|) and vary 
both Q 2 and M 2 . The function C(Q 2 ) (and consequently 
gxif>L>(Q 2 )) should not depend on M 2 , so we limit our fit 
region to 3.0 GeV 2 < M 2 < 3.5 GeV 2 where C{Q 2 ) is 
clearly stable in M 2 for all values of Q 2 . 



8ir 2 m 2 x 



Blj — >TT7T7tIuJi ("72) 



We do the fitting for s 



1/2 



4.4 GeV as the results 



do not depend much on this parameter, the results are 



which, together with Eq. (fTTj) gives us: 

r (X -> J/tPtt + tt-tt°) = (9.3 ±6.9) MeV. (73) 

The result in Eq. (|73[) is in complete agreement with 
the experimental upper limit. It is important to notice 
that the width grows with the mixing angle 9, as can be 
seen from Eq. (|68p . while the mass grows with 9. There- 
fore, there is only a small range for the values of this an- 
gle that can provide simultaneously good agreement with 
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the experimental values of the mass and the decay width, 
and this range is 5° < 9 < 13°. This means that the 
X(3872) is basically a cc state with a small, but funda- 
mental, admixture of molecular DD* states. By molecu- 
lar states we mean an admixture between D°D*°, D°D*° 
and D+D*-, D~ D*+ states, as given by Eq. (J^J). 

V. CONCLUSIONS 

We have presented a QCDSR analysis of the two-point 
and three-point functions of the A"(3872) meson, by con- 
sidering a mixed charmonium- molecular current. We find 
that the sum rules results in Eqs. (f35|) and (f73|) are com- 
patible with experimental data. These results were ob- 
tained by considering the mixing angle in Eq. (j 10)) in the 
range 5° < 9 < 13°. 



We have also studied the mixing between the 
D°D*°, D°D*° and D+D*-, D- D*+ states by impos- 
ing the ratio in Eq. ([1]). In accordance with the findings 
in ref. we found that the mixing angle in Eq. (|59[) is 
a ~ 20°. 

With the knowledge of these two mixing angles we con- 
clude that the X(3872) is basically a cc state (~97%) 
with a small, but fundamental, admixture of molecular 
D°D*°, D°D*° (-88%) and D+D*-, D~ D*+ (-12%) 
states. 

This small molecular component could, in principle, 
be a consequence of neglecting the two-hadron reducible 
contribution in the phenomenological side. However, as 
argued in section III, we expect the 2HR contribution to 
be small and the results to hold even if we had taken it 
into consideration. 
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